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Abstract: Using generalized Konishi anomaly equations, it is known that one can ex- 
press, in a large class of super symmetric gauge theories, all the chiral operators expecta- 
tion values in terms of a finite number of a priori arbitrary constants. We show that these 
constants are fully determined by the requirement of gauge invariance and an additional 
anomaly equation. The constraints so obtained turn out to be equivalent to the extremiza- 
tion of the Dijkgraaf-Vafa quantum glueball superpotential, with all terms (including the 
Veneziano-Yankielowicz part) unambiguously fixed. As an application, we fill non-trivial 
gaps in existing derivations of the mass gap and confinement properties in super Yang- Mills 
theories. 
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1. General presentation 



The aim of this paper is to present a general proof of the Dijkgraaf-Vafa quantum equations 
of motion [1], 

that determine the glueball superfields expectation values sj = (Sj) in a large class of 
J\f = 1 super symmetric gauge theories. The holomorphic function Wy)y{si) is the sum of 
a suitable Veneziano-Yankielowicz term [2] and of an infinite power series in the sjs. The 
equations (1.1) are generally interpreted in a low energy effective action framework: the 
gauge invariant composite fields Si are assumed to be good low energy degrees of freedom, 
with low energy effective superpotential given by Wdv- The conditions (1.1) are then the 
dynamical equations of motion in the low energy effective theory. This interpretation is 
extremely difficult to justify from first principles in the gauge theory, as stressed for example 
in [3]. Our main result is to show that (1.1) actually follows from general consistency 
conditions in the gauge theory, in some sense analogous to the Dirac quantization condition. 
More precisely, we shall demonstrate that the vanishing of the differences {dsj — c?s^)Wdv 
follows from gauge invariance, and that the vanishing of the sum sjOsjWdy follows from 
a Ward identity. We thus interpret the equations (1.1) as being non- dynamical identities. 

In the dual closed string theory formulation [4] , Wby is a genuine effective superpoten- 
tial associated with the non-trivial background fluxes [5]. Prom this point of view, deriving 
(1.1) amounts to proving that the closed string background proposed in [4], and their gen- 
eralizations based on the geometric transition picture (see [6] for non-trivial examples), are 
indeed the correct duals of the gauge theory. We thus have a dual interpretation of the 
equations (1.1): they are non-dynamical from the gauge theory, or open string point of 
view, and dynamical from the closed string point of view. This is a nice new example of 
a familiar phenomenon for dual formulations of the same theory (the basic example being 
the electric-magnetic duality of the Maxwell's equations, under which the Bianchi identities 
and the dynamical equations of motion arc exchanged). 

A startling feature of the equations (1.1) is that they provide a unifying framework 
to derive essentially all the known exact results in supersymmetric gauge theories. For 
example, they have been used to obtain the Seiberg-Witten solutions of A/" = 2 theories [7] , 
to study Seiberg dualities [8] , or to find new results on the quantum space of parameters in 
M =1 theories [9, 10, 11, 12]. In some cases, the same results can also be obtained by using 
various techniques (S-duality, singularity theory, integrable systems, geometric engineering, 
mirror symmetry, brane constructions, geometric transitions etc. . . ), but at the expense 
of making several non-trivial assumptions. Actually, there are very few derivations from 
first principles. The one exception is Nekrasov's explicit summing-up of instanton series 
[13], which yields the Seiberg-Witten prepotentials in A/" = 2 super Yang-Mills theories. 
The main application of the proof of (1.1) is to provide automatically assumption- free 
derivations from first principles of all the above mentioned exact results. For example, the 
derivation of the Seiberg-Witten N = 2 prepotentials along these lines do not even assume 
that only instantons contribute, as in Nekrasov's, put prove it (see Section 5). Another 
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important application is to the mass gap and confinement problems in = 1 super Yang- 
Mills. The main argument in favor of confinement was given by Seiberg and Witten in 
their classic papers [14]. The idea is to perturb a J\f = 2 theory by adding a tree level 
superpotential. The argument relies on two main hypotheses: the solution of the M = 2 
theory is correct; the M = 1 theory obtained by perturbing M = 2 creates a mass gap. We 
are able to waive these hypotheses in Section 5. 

The paper is organized as follows. In Section 2, we review some useful concepts and 
results, including the work of Cachazo, Douglas, Seiberg and Witten [3] on generalized 
Konishi anomaly equations. These equations are used in Section 3 to study, for any adjoint 
field X in the theory, the quantum characteristic function 

^{z) = (det(z - X)) . (1.2) 

Our strategy is to show that gauge invariance puts general constraints on ^ that turn 
out to be equivalent to {dgj — dsj)W-DY = 0. In Section 4, we find a Ward identity that 
provides the missing equation in the form SidsiWuY = 0. We also comment on the 
possible existence of Kovner-Shifman vacua [15]. Finally in Section 5 we apply the results 
to the mass gap and confinement problems, giving en passant the derivation of the Seiberg- 
Witten solution for J\f = 2. Our arguments are completely general, but for concreteness we 
mainly focus on the M = I theory with gauge group U(Ar), one adjoint matter superfield 
X and ATf < 2N fundamental flavors. We also give details on the case with two adjoint 
fields, because it illustrates nicely some non-trivial aspects of the proof. 



2. Introductory material and discussion 
2.1 The quantum effective superpotential 

A most important object in four dimensional N = 1 supersymmetric gauge theories is the 
quantum effective superpotential M^iowj which determines the F-term part of the quantum 
effective action. It is defined by a path integral over the vector and ghosts supermultiplets 
V and the matter chiral multiplets X, with boundary conditions determined by the choice 
of a supersymmetric vacuum |0), 

^i/d4x(2iVRe/d2ewl°^(rK)+D-terms) ^ f ^^^^i ^ A-^x J^{V,X;Tk) _ (2.1) 

J\Q) 

The normalization factor N is chosen for convenience to coincide with the number of colors 
of the gauge group U(Ar). The Tk are external chiral superfields, with lowest components 
tx, that couple to gauge invariant operators Ok in such a way that 

(0|O^|0) = ^. (2.2) 

In general, the knowledge of Wiow is equivalent to the knowledge of the full set of chiral 
operator vacuum expectation values in the theory. 
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For example, consider the U(A^) gauge theory with one adjoint X, N{ < 2N flavors of 
fundamentals Qf and anti-fundamentals Q-^ chiral multiplets, and tree level lagrangian 



^= -^Im 
47r 



niv I d^9 TrpW^Wc 

+ 2NRe J e (Ttf WtreeiX) + m /' {X)Qf,^ +D-teims. (2.3) 



The traces Trp are in the fundamental representation. The bare gauge coupling ruv at 
scale Auv is given in terms of the dynamically generated scale A by 

/ A \2N-N( 

2i7rruv = ^'^(x^J^J ■ (2-4) 

The one exception is when A'f = 2N, in which case ruv = t is not renormalized. The 
derivative of the tree-level superpotential is taken to be of the form 

d d 

<ee(^) = Y.9kX^ = J{{X - ai) . (2.5) 

fe=0 /=! 

The polynomials 

mf{X) = ^rn/gX" (2.6) 

q>0 

are such that 

e 

det f,f'm{X) = Uol{{X-bQ). (2.7) 

Q=l 

The most general vacuum 

\Ni,ki;...;Nd,kd;i^i;...;iye) (2.8) 

is specified by the numbers Nj and uq, uq = or I, of eigenvalues of the matrix X that 
classically are equal to a/ and bq, with pattern of gauge symmetry breaking U{N) 
'[J{Ni) X • • • X \J{Nd), and by the integers kj, defined modulo ZjVj, that correspond to low 
energy chiral symmetry breaking. In particular, 

d e 

^Ni + J2^Q = N. (2.9) 

1=1 Q=l 

In any vacuum \Nj,ki\VQ), the relations (2.2) read 

-l-(TrpX'=+i> = ^, fQ^X«Q/'> = ^> (2.10) 
k + l^ ' dgk ^ ' ' dm 

1,1 



where we have introduced the glueball superfield 



(2.12) 
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When N{ = 2N, InA^^"^' is replaced by 2i7rr in (2.11). 

By a standard non-renormalization theorem, VFiow is not corrected in perturbation 
theory, but there are in many cases non-trivial non-perturbative corrections, typically given 
by infinite sums over (fractional) instantons. The calculation of Wiow is thus a difficult 
strong coupling problem. 

2.2 The Dijkgraaf-Vafa approach 

Dijkgraaf and Vafa made the remarkable conjecture [1] that the calculation of Wiow iii ^ 
large class of A/" = 1 supersymmetric gauge theories could be reduced to a simple recipe. 
The recipe goes in two separate steps. 

2.2.1 Step one: the matrix model 

First one introduces a zcro-dimensional gauged matrix model whose potential coincides 
with the tree-level superpotential of the gauge theory. For example, for the model (2.3), 
we use the one-matrix model coupled to flavors, for which the partition function is an 
integral over a n x n hermitian matrix X and n- vectors and Q f , 

^n-r/s- = ^1 dXdQdge-?^^-^*-(^)-?"'5^-/'(^)Q/' . (2.13) 

The normalization factor ^ is independent on the couplings 

tK = {gk,m/j (2.14) 

in the tree-level action and will be discussed later. The size n of the matrix and the vectors 
is a dummy variable and is not related to the number of colors N in the gauge theory. We 
are actually interested in the n — > oo limit, or more precisely in the sum over sphere and 
disk diagrams around a particular classical (s 0) solution, which yields an expansion of 
the form 

= J^s + ^^d + C»(l/n2) . (2.15) 

The classical solutions are classified by the integers nj and z/g, = or 1, where n/ and 
VQ give the number of eigenvalues of X equal to aj and bq in the s — > limit. In the 
n — > GO limit, the partition function, obtained by expanding around a particular matrix 
model vacuum 

(2.16) 

is a natural function of the filling fractions 

si = s— (2.17) 
n 

and of the vq. Actually, since the sphere partition function J^s is insensitive to the coupling 
to the fundamental flavors, it depends on the sj only, whereas the disk partition function 
J-'^ depends on both the sj and the i^q. 

We then introduce the Dijkgraaf-Vafa superpotential [1, 16] 
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The "weak" Dijkgraaf-Vafa conjecture states that 

where the vq parametrizing the matrix model (2.16) and gauge theory (2.8) vacua are taken 
to be the same. The equations (2.19) determine an infinite set of chiral operators vacuum 
expectation values in terms of a finite number of undetermined constants sj. These matrix 
model filling fractions are identified in the gauge theory with the gluino condensates in the 
V{Ni) factors of the low energy gauge group, 

SI = {Si) , (2.20) 

or are taken to be zero if Ni = (as explained in [3], the operators Si are defined as in 
(2.12) with the insertion of a suitable projection operator). 

2.2.2 Step two: the extremization of VFdv 
The "strong" Dijkgraaf-Vafa conjecture states that 

W^iow = WMs^i) , (2.21) 

where the are the particular values of the si that solve the quantum equations of motion 
(1.1). 

The function Wby splits into a part Wyy that contains all the ii^-independent terms 
in VFdv and the tii--depcndent part that we denote PFmm, 

Wovisi; tx) = Wyyisi) + Wmm{si; Ik) , (2.22) 

It is important to realize that the weak Dijkgraaf-Vafa conjecture does not depend on the 
form of Wyy, since any fft-independent term drops from (2.19). On the other hand, Wyy 
plays a crucial role in the equations (1.1). A central result of the present work will be to 
derive, in Sections 3 and 4, the form of Wyy from first principles. 

With our definition (2.13), l^vY depends on the normalization factor A priori, this 
normalization factor may be a complicated function of s and of the d3mamically generated 
scale A. For example, for the model (2.3), we may postulate ^ oc A^s** with a + 36 = 
— — 2nArf for dimensional reasons. The correct choice (that we will be able to justify 
later) turns out to be 

^ oc Z^-nHnNtg-nNt _ ^2.23) 

There can be an additional purely numerical factor (independent of the tx, si or A), that 
can be used to set a convenient normalization of A. 

A last important comment is that with the choice (2.23), (2.13) and (2.18) would imply 
that Wyy = — '^j Nisilnsi, which is slightly misleading. The correct formula is 

WVy(s7) = -5I*-flnsf^ (2.24) 
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With this prescription, solving (1.1) automatically yields the full set of chirally asymmetric 
vacua (2.8) labeled by the integers kj. Let us note that Wyy is the sum of the independent 
pure gauge theory contributions postulated in [2] for each factor \J{Nj) of the low energy 
gauge group. For this reason, it is often called the Veneziano-Yankielowicz part of the 
Dijkgraaf-Vafa superpotential. 

2.3 The generalized Konishi anomaly equations 

The weak form of the Dijkgraaf-Vafa conjecture was explained in a seminal paper by Cac- 
hazo, Douglas, Seiberg and Witten [3, 17].^ These authors showed that the loop equations 
of the matrix model, that determine all the matrix model correlators as a function of 
the filling fractions si, were in one-to-one correspondence with a set of generalized Kon- 
ishi anomaly equations in the gauge theory, the precise form of the correspondence being 
equivalent to (2.19). 




Figure 1: The non-compact two-sheeted Riemann surface C, with the contours ai, (3i and 777 

used in the main text. The contours /?/ are defined modulo an integer multiple of a/. The points 
2; = /Lto on the first and second sheets are taken to infinity. 

The consequences of the generalized Konishi anomaly equations, and in particular of 
the equations (2.19), for the model (2.3) can be summarized as follows. Let us introduce 
the generating functions for the (TVf^''), (Trp W"WaX^) and (^Q^'X^Qf), 

1 



y{z) 

and the convenient combination 



Tr 



z-X 
1 / 



:— (Trp " 

167r2Ar\ z-X 



z-X 



Qf 



n^) = Wl,^{z)-2S^{z). 



The function W satisfies a quadratic equation 



A(^), 



(2.25) 
(2.26) 
(2.27) 

(2.28) 
(2.29) 



^Unfortunately, a proper non-perturbative discussion of the arguments in [3] has not yet appeared in 

the hterature. In the present paper we focus on the strong form of the Dijkgraaf-Vafa conjecture, assuming 
that the weak form of the conjecture, or equivalently the results of [3], are correct. 
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where A is a polynomial of degree deg A < d — 1 = dcgW^j.^^ — 1. In a vacuum with low 
energy gauge group U(A^i) x • • • x U(A''^), 1 < d < d, the following factorization in terms 
of a degree d — d polynomial H^_^ takes place, 

d 

Wll^iz) - A{z) = Hl_lz) \{{z- a7)(z - a]) , (2.30) 

7=1 

and thus ^ and are meromorphic functions on the genus d — 1 hyperelliptic surface 

d 

C: y^ = l[{z-aj){z-ai). (2.31) 
1=1 

This two-sheeted Riemann surface, with some contours defined on it, is depicted in Figure 
1. The polynomial A, and thus the surface C, is completely determined by the factorization 
condition (2.30) up to d parameters that can be conveniently chosen to be the aj-period 
integrals of S^dz, 

i^£^d. = »,. (2.32) 

In particular, we have 

Y^si = s = {S), (2.33) 

the expectation value of the glueball supcrficld defined in (2.12). In the classical limit, 
sj for all / and we have aj ~ a^^ ~ a/. Equations (2.28) and (2.29) also show that 
=5^(z) is the sphere matrix model resolvent, 

y{z) = ((-Tr^^—} , (2.34) 

W n Z — X // sphere 

where the double bracket notation is used for the matrix model expectation values (which 
are clearly different in general from the expectation values in the gauge theory). The 
parameters sj introduced in (2.32) are thus the same as the filling fractions (2.17) of the 
matrix model. 

f' 

The generating functions ^^^-^ and ^ are also meromorphic functions on the curve 
(2.31). We shall use the gauge theory resolvent ^ extensively in the following, so let us 
describe its structure in details. Let us introduce 

Mz) = ■ (2.35) 

It is straightforward to show that the il)idz form a basis of log-normalizable holomorphic 
one-forms on C, 

V-/ = ^ • (2.36) 

y 

The polynomials pi = z'^~^ + • • • arc fixed by the constraints 

/ ^Pidz = dij. (2.37) 
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The gauge theory resolvent is then given by 

M{z) = NiM^) + gd{z) , (2.38) 

where is the disk matrix model resolvent (obtained by summing disk diagrams only), 

If we define the "classical" sheet of the Riemann surface (2.31) by the condition 

AT 

^(z) ~ — , (2.40) 

z^co Z 

then M has poles on the classical sheet at z = bq with residue uq li uq = \, and poles on 
the second sheet at 2; = 6q with residue 1 — if i/q = 0. We also have the useful identities 

— / ^{z) dz = Ni. (2.41) 

For completeness, let us finally mention that 

The above equations determine all the (TrpX'^), {WWaTrpX'') and {^Qf'XiQf) in 
terms of the sj. Note that taking into account (2.35), (2.38) and (2.42) are perfectly 
consistent with (2.19) and (2.18). 

2.4 Another formulation of the quantum equations of motion 

A last important technical result that we need is an elegant reformulation of the quantum 
equations of motion (1.1) as constraints on the period integrals of ^dz [12]. Using the open 
contours /?/ defined in Figure 1, and with the understanding that /xq is taken to infinity, 
the equations (1.1) are equivalent to 

1 /■ . , 1 . f7o_A'^-^' 



2zw.^^^^^^ = 2^;^^^T^ + ^^' ^'-^^^ 



where the kis are integers defined modulo Ni. In particular, the periods on the compact 
cycles 7/j = (3i — (3j are integers, 

I M{z) dz = ki-kjeZ. (2.44) 

Since this result is crucial for our purposes, we now present a short proof. The starting 
point is the so-called special geometry relations of the matrix model, 

^'•^^ ^ Vjdx + 21n^, (2.45) 



dsidsj Jfij A 
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as well as similar equations for the derivative of the disk partition function, 




dx + ln(//Q 



/Uo). 



(2.46) 



dsi 



The details on the derivation of these matrix model relations can be found for example in 
[18]. Equations (2.18) and (2.38) then yield 



This would correspond to a Veneziano-Yankielowicz term —Ni Ins/ in dWnY / dsj, whereas 
the correct term is — Ins^^ as in (2.24). Taking into account this subtlety, we obtain (2.43). 
The ki are defined modulo Ni because the contours /?/ are defined modulo an integer 
multiple of a/ (the transformation /?/ /?/ + a/ is induced by continuously permuting aj 
and a\ , which obviously does not change the curve C nor any gauge theory correlator, and 
thus yields the same solution). 

The most natural form of (2.43) is 



Prom the above discussion, we know that there are automatically n7=i distinct solutions 
to these equations, labelled by the ZjVj-valued integers kj. 

2.5 Discussion 

The Konishi anomaly equations can be generalized to a very large class of models and 
put very powerful constraints on the set of chiral operators expectation values. However, 
any choice of the parameters si yields a perfectly sensible solution of the equations. For 
instance, the choice = for all I corresponds to no quantum corrections at all to VFiow 
Clearly, the anomaly equations do not address directly the problem of the strongly coupled 
gauge dynamics. This was explained in details already in the original paper [3]. Neverthe- 
less, these equations do provide some important indirect information, since we know that 
the true answer must correspond to some particular values of the s/. In particular, we 
know that and ^ are functions defined on the two-sheeted surface (2.31), since this is 
true for all values of the sj. The knowledge of this analytic structure will be of great help 
in the following. 

Unlike the Konishi anomaly equations, the strong part of the Dijkgraaf-Vafa conjecture 
deals directly with genuine non-perturbative effects in the gauge theory, including infinite 
series of (fractional) instanton contributions and chiral symmetry breaking. It is also closely 
related to the problems of the mass gap and confinement in A/" = 1 theories, as we shall 
review in Section 5. 

The original formulation of the conjecture in terms of the extremization of the function 
Wdv suggests that Wdv is a low energy cfi'cctivc supcrpotential and the equations (1.1) 
purely dynamical. It would clearly be very difficult to devise a proof along these lines. 




(2.47) 




(2.48) 
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For example, there is no convincing argument showing that the glueball fields Sj are the 
good low energy variables. Actually, this interpretation is logically disconnected from the 
Dijkgraaf-Vafa conjecture, and it is even unclear whether it is fully correct or not. 

The equations (2.43) or (2.48), even though they are mathematically strictly equivalent 
to the original Dijkgraaf-Vafa formulation in terms of a superpotential, are beautifully 
simple and offer an attractive alternative starting point. At the very least, they suggest 
an entirely different interpretation. In particular, the equations (2.44) look like simple 
quantization conditions, which are often derived from general consistency constraints. We 
show in the next Section that this is indeed the case, (2.44) being a consequence of gauge 
invariance. The full set of equations (2.48) then follows from an additional constraint, 
which is a simple Ward identity that we derive in Section 4. 

The relation between (2.44) and gauge invariance is not totally unexpected. Assuming 
that confinement is realized in the U(iV/) factors of the low energy gauge group through 
the condensation of dyons of unit magnetic charge (the usual 't Hooft's mechanism), one 
can show that the integers kj in (2.43) coincide with the electric charges of the dyons [11]. 
The Dirac-Schwinger-Zwanziger quantization condition, which is well-known to be a direct 
consequence of gauge invariance, is then equivalent to (2.44). Our main task in Section 3 
is to present a simple derivation of the same result that is free from any assumption on the 
low energy behaviour of the theory. 

2.6 Other models 

An important property of (2.44) is its universality. Even though the details of the Dijkgraaf- 
Vafa prescription (matrix model, definition of Wdv) are model-dependent, the extremiza- 
tion of the Dijkgraaf- Vafa superpotential always yields quantization conditions of the form, 
(2.44), with ^ being defined as in (2.25) to be the resolvent associated with any adjoint 
chiral superfield in the theory. This fact will find a simple explanation in the next Section, 
but we wish to illustrate it using elementary methods on a representative example. 

Let us introduce the theory with two adjoints X and Y and a tree-level superpotential 
of the form 



This case is particularly interesting, because there are now two natural gauge theory resol- 
vents. 



for which constraints of the form (2.48) can be expected to be valid. This is non-trivial, 
because the constraints are then redundant, and not obviously consistent with each other. 
The matrix model we have to consider is 



Wuee{X, Y) = ViiX) + V2iY) - nXY , V( = P,V^ = Q. 



(2.49) 




(2.50) 




(2.51) 



with a Dijkgraaf-Vafa superpotential 




(2.52) 
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It is then possible to show that, provided the normahzation constant is chosen to be 

jY oc A-"V25-nV2 ^ (2.53) 

the equations of motion (1.1) take a form similar to (2.48). In particular, the quantization 
conditions ^ 

— I ^^{z)dz e Z (2.54) 

and 

i M^{z)dz G Z (2.55) 

turn out to be equivalent to each other. The choice (2.53) for jY is also necessary for the 
standard interpretation in terms of Veneziano- Yankielowick terms to be valid. The contours 
^fj and ^Jj are similar to the contours 77 j of Figure 1, but with the 7^ (respectively '^Jj) 
going through the cuts of the resolvent ^^-^ {^^)- The proof of (2.54) and (2.55) uses the 
generalized Konishi anomaly equations (they have not appeared in the literature for the 
model (2.49), but they are straightforward to derive) and the special geometry relations 
for the matrix model that are derived in [6]. Full details are given in [18]. 

3. The quantization conditions 
3.1 Presentation of the problem 

Let M be an arbitrary adjoint chiral superfield. For example, in the model (2.3), we may 
consider any power of the elementary field X\ in the model (2.49), any combination of the 
form Y\.k XP''YP'<= would qualify. The resolvent for M is defined by 

The classical resolvent has the form 

^cl(^)=E^^' M 

where the integers Nj, -/V/ = N, give the number of eigenvalues of M that are equal 
classically to aj. Quantum mechanically, we assume that ^ is a multi-sheeted analytic 
function. On the classical sheet, which is defined by the property that the large z asymp- 
totics is the same as in the classical theory, 

N 

^(z) ~ — , (3.3) 

z— »oo z 

some of the poles (say the a/ for 1 < / < d) in (3.2) are replaced by square root branch 
cuts, and the others (the aq fov Q > d + I) remain poles that may be shifted from their 
classical value to a quantum value 6q. The analytic structure is thus similar to the one 
given in Figure 1 (in general there could be additional branch cuts on the second sheet, and 
the classical sheet could be glued to many different sheets by the square root branch cuts; 
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we shall work with an analytic structure of the form given by Figure 1, because the other 
cases amount to straightforward generalizations). As explained in Section 2, the above 
assumptions on ^ can be derived from the generalized Konishi anomalies in many cases 
(including the cases that are relevant for our purposes), and so they are not restrictive. 
The classical limit (3.2) and the analytic structure of Si imply trivially that 



2^7^ 



z) d^ = A^/ G Z . (3.4) 

Otl 

We want to prove that gauge invariance implies that the other period integrals of MAz 
must be quantized as well, 

^ j, ^{z)dz G Z. (3.5) 

3.2 On gauge invariance 

In [3], the single-trace gauge invariant chiral operators in the theory with no flavors were 
classified (a similar classification can be achieved in many other models). The full set is 
given by TrpX^, Tip X'^W'^ and Tip X'^W'^Wa- We call the ring generated by these ele- 
ments the reduced chiral ring of the theory. Because of the factorization of chiral operators 
correlation functions, the structure of this ring is fully determined by the knowledge of the 
functions ^{z) and y{z) defined in (2.25) and (2.26). 

The most general element of the reduced chiral ring can be written as a finite sum 
of finite products of single trace chiral operators. This does not exhaust the full set of 
interesting chiral operators. In particular, ,'^(z) and '^(z) themselves are not in this ring. 
More generally, it is natural to consider operators of the form 

0{z) = J2 Okz-" , (3.6) 

k>K 

where the Ok are elements of the reduced chiral ring. The series defining the expectation 
values 

ff{z) = ^ (3.7) 

k>K 

are assumed to have a non-zero radius of convergence around z = oo. Operators like 
(3.6) are usually introduced because they are generating functions for the (Ofe). This 
interpretation comes from the expansion at large z. In particular, when z is in the disk 
on convergence, there is little to learn in considering 0{z) instead of the individuals O^. 
However, we shall see that considering a priori 0{z) for any value of z introduces some 
interesting subtleties. 

When z is not in the disk of convergence, i^{z) is defined by analytic continuation. 
The result is in general a multivalued function on the z-plane. This multivaluedness is not 
obviously consistent with gauge invariance. Gauge invariance is manifest in an expansion 
like (3.6), but this expansion is only valid on a single sheet, and in general a similar 
expansion does not exist on the other sheets. To understand how the problem can arise, 
let us focus on the case where ff{z) is built entirely in terms of a single adjoint field M in 
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the theory. The operators are then polynomials in the Trp Af, and thus ^{z) depends 
only on the eigenvalues rrii of M on the sheet where (3.7) is valid, and thus on all the other 
sheets by analytic continuation, 

ff{z) = ff{z\mi,...,mN). (3.8) 

The fundamental constraint from gauge invariance is then that ^{z;mi) must be a sym- 
metric function of the mi. If a- is the operator implementing the gauge transformation 
associated with the permutation a G Sn C U(A'^) of the eigenvalues, we must have 

{a ■ ff){z; mi) = ffiz; m^(i)) = €?(z; m^) . (3.9) 

This constraint is obviously satisfied when the series (3.7) converges, but is not insured for 
the analytic continuations. 

To illustrate this point, let us consider an extremely simple classical toy example, 

0{z) = (ttf Vz - mY = N^z - N Trp M + ^{{TiF - N TiF M^) + ■ ■ ■ . (3.10) 

In terms of the eigenvalues, we have 

^{z; mj) = ^ ^ {z - mi){z - ruj) . (3.11) 

The function i^{z;mi) has square root branch cuts that join pairs of eigenvalues. Let us 
start from a value of z for which the series in the right hand side of (3.10) converges, 
and let us consider the analytic continuation through the cut between and nij^^ . This 
analytic continuation amounts to changing the sign of the term \/{z — miQ){z — mj^) in 
(3.11), and is thus clearly inconsistent with the symmetry of & under the permutation of 
rriiQ (or m^g) with any other eigenvalue when i {io) jo}- This shows that 0{z) defined 
by (3.10) is not a good gauge invariant operator for all values of z. We could of course 
compute it at large z in the sheet where it is gauge invariant, and use it as a convenient 
book-keeping device for the Trp X^. However, considering (Trp Vz — Af) for all values of 
z is inconsistent with gauge invariance. 

In the general case, let us label by an index a > the different sheets of the analytic 
function ff, and write Ga{z\mi) its value on the a**^ sheet. Suppose that the 0**^ sheet 
contains the disk of convergence of (3.7), and thus in particular that 

{g ■ ^o)(^; mi) = ffo{z; rrii) . (3.12) 

If the action of the operator a- is smooth, the same equation must automatically be valid 
on any other sheet, by analytic continuation. This shows that if there exists a sheet, say 
number a, for which the permutation symmetry is violated, 

{a ■ ^a){z;mi) ^a{z;mi) , (3.13) 

then the gauge transformation a- must be singular when acting on ^{z) for some values 
of z. For example, if the 0*'^ and a^^ sheets are joined by a branch cut, then a- is clearly 
not well defined (two- valued) on this cut. This inconsistency is similar to the problems 
that appear in the presence of a Dirac string when the Dirac quantization condition is not 
satisfied. 
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3.3 The quantum characteristic function 
3.3.1 Generalities 

Let us introduce the characteristic function of M, 



Fi = Trp M , 
^2 = ^ Trp - 

Picking a particular value of N 



^{z) = (det(z - M)> , 



(3.14) 



with the basic property that 

5[^=«W. (3.15) 

By using 

det(z - M) = z^e^^ H^-M/z) ^ ^^^^-^ 
we can obtain an expansion of the form (3.6), 

det(z - M) = - ^ (3.17) 

where the elements of the rcdTiccd chiral ring Fk = I Trp + ■■■ are homogeneous 
polynomials in the Trp M'^ of degree k (M being of degree one) whose form does not 
depend explicitly on the number of colors N . For example, 



-1(TVFM)^ (3.18) 
- ^ (Trp M) (Trp M^) + i (TVp M)\ etc. . . 
amounts to imposing relations in the chiral ring, 



Fjv+g = 0, q>l, (3.19) 

that arc equivalent to the fact that classically det(z — M) is a polynomial of degree N . 
Using (3.19) recursively, one obtains the standard indentities 

TVp M« = P,,ci (TVp M, . . . , IVp M^) , q>\, (3.20) 

where the ci ^ire homogeneous polynomials of degree q depending only on the N inde- 
pendent Casimh Trp M'', l<k<N, of U(iV). 

The chiral ring relations (3.19) or (3.20) can get non-perturbative corrections, 

TVpM'' = Pg(TrpM,...,TrpM^;A) . (3.21) 

The generalized Konishi anomaly equations put strong constraints on the form of the 
polynomials Pq. In particular, if the filling fractions sj are not all equal to zero, the 
quantum corrections must be non-trivial. The equations (3.21) then imply that the quan- 
tum characteristic function (3.14) or (3.17) is not a polynomial, but an infinite series, the 
non-polynomial terms precisely generating the quantum corrections to the chiral ring. 
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3.3.2 The integral representation 

The characteristic function ^(z) is thus multivalued (this follows from (3.15) and the 
multivaluedness of ^{z)). The fact that a correlator can be multivalued in a field theory is 
of course not an inconsistency, and actually is a very common phenomenon. For example, it 
happens when the theory has different branches, or loosely speaking different phases. More 
concretely, when we consider the path integral for a theory with a (discrete or continuous) 
moduli space of vacua, this path integral is inherently multivalued, the ambiguity coming 
from the ambiguity in the choice of vacuum. It is also well-know that, in some cases, 
the different branches can be continuously related to each other by analytic continuation 
[9, 10, 11, 12]. We are going to explain that the multivaluedness of ^{z) can be understood 
in similar terms. Note however that there is an important difference: in cases like ^{z), 
for which z enters the path integral only through the insertion of a local operator, the 
boundary conditions at infinity cannot depend on z. This means that the gauge theory 
vacuum (2.8) must be the same for all branches of 

The characteristic function can be expressed, in a manifesly gauge invariant way for 
all values of integral over anticommuting variables, 

det(^ - M) = y d^dx exp [xa [6^ - z-'M%) . (3.22) 

Strictly speaking, this formula is only valid classically. Quantum mechanically, a smearing 
in space-time is necessary, and the integral in the right hand side of (3.22) is to be inter- 
preted as being a path integral localized around a certain space-time point x, in the limit 
where the smearing distribution goes to the Dirac 5. This smearing is very similar in nature 
to the smearing of the glueball field introduced in [19]. In particular, (xaM'^j^ip^)^ can be 
non-zero for arbitrary powers k. When expanded around ip^ = Xa = 0, the integral (3.22) 
then yields an infinite series in z. Using Wick's theorem and the contraction XaV^ = 
the limit where the smearing is removed, it goes to the expansion (3.17). 

Expanding aroimd = Xa = is not the only possibility. Indeed, the classical 
equations of motion for (3.22) read 

Ml^'' = zr , XaM% = zxb . (3.23) 

When z coincides with an eigenvalue of M, there is a moduli space of solutions ("Higgs 
branch") parametrized by the gauge invariant "meson" 

a = Xar ■ (3.24) 

This is a hint that (3.22) can indeed be multivalued. Quantum mechanically, the eigenvalues 
of M fiuctuate. On the one hand, these fiuctuations can lift the fiat direction in a, and on 
the other hand the "Higgs branch" can exist for any z, with a given value of {a) = aii{z). 
This behaviour is reproduced for example by a simple toy model for which the quantum 
corrections are modeled by a gaussian integration over M (note however that the gauge 
theory fluctuations arc much more subtle: single trace operators behave as if the eigenvalues 
were fluctuating, whereas multi-trace correlators factorize as if they were dominated by a 
non- fluctuating master field). 
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Taking the derivative with respect to z, the formula (3.22) shows that the condensate 
(a) is nothing but the gauge theory resolvent 

{a){z)=^{z). (3.25) 

We thus have a nice interpretation of the multivaluedness of ^{z), that correspond to the 
different "phases" of the integral (3.22). 

3.3.3 A simple proof of the quantization conditions 

The above interpretation may be innocent-looking, but it has deep consequences. The path 
integral (3.22) can be multivalued because it can have several "vacua." In particular, once 
we pick a vacuum, all the ambiguity must be lifted. But the vacua arc characterized by 
the set of gauge invariant expectation values that can be built from the elementary fields 
over which we integrate. In the case of (3.22), all we have is the gauge invariant meson 
(3.24), obtained by taking the derivative with respect to z, or higher powers of a, obtained 
by taking higher derivatives. This implies that, even though (3.22) may be a multivalued 
function of z, it must have precisely the same sheets and branch cuts as {cr){z) = ^(z). In 
other words, it must be a meromorphic function on the Riemann surface (2.31). 

Let us analyze the consequences of this fact, focusing on the case of the theory (2.3) for 
concreteness. By integrating the relation (3.15) from the point at infinity on the classical 
sheet to z, we obtain 

^{z) = exp f ^{z') dz' . (3.26) 

This formula shows that ^ is generically a multivalued function on the curve (2.31), 
because one must specify the contour from the point at infinity fXQ to z to do the integral. 
A first simplification occurs thanks to the conditions (3.4). We can define unambiguously 

^{z) on the classical sheet by integrating along any path that remains on this sheet. The 
answer does not depend on how many times the path circles around the branch cuts of 
This means that ^ and ^ have the same branch cuts on the classical sheet, and the latter 
can be equally well characterized by the properties (3.3) or by 

^(z) ~ z^ . (3.27) 

Z—KX> 

We can say more about the classical sheet branch cuts of as for they are square 
root branch cuts. This is so because going twice across the cut corresponds to going along 
a contour whose homotopy class is trivial. Let us also note that the poles of ^ become 
simple zeros of Complications arise when we consider the analytic continuations across 
the classical sheet branch cuts. Unlike for if wc go through two different branch cuts, 
we obtain a priori two inequivalent continuations of Equation (3.26) shows that this 
does not occur if and only if 

exp * M{z) dz = l, (3.28) 

in which case =^ and ^ are defined on the same surface (2.31) as required. Obviously, 
(3.28) is equivalent to the quantization conditions (3.5). 

We are now going to show that the same result also follows from a rigorous implemen- 
tation of the constraints (3.9) on 
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3.4 The permutation symmetry 
3.4.1 Dependence of the eigenvalues 

Wc now focus on the model (2.3) with M = X for concrctcness. Let us go to a gauge where 
X is diagonal (this is the supersymmetric version of 't Hooft's abelian projection [20]), 

X = diag(a;i ,...,xn)- (3.29) 

The gauge fixing procedure in the gauge theory, which mimics tlic gauge fixing in the 
associated matrix model (2.13), was discussed in [21]. The eigenvalues xj are operators 
in the quantum theory. In a given vacuum (2.8), we expand each eigenvalue around its 
classical value. For example, we may choose 

Xi = ai + Xi ioT ^Nj + 1 <i <^Nj . (3.30) 

J<7 J<I 

The gauge theory resolvent can be decomposed in a natural way as 

d e 
M{z) = ^ Ni^iiz) + (3.31) 

1=1 Q=l ^ °^ 

with 

1 / dz' !M{z' 



E/. + l<.<EiV.. (3.32) 

The variable z in (3.32) is taken to be outside the contour aj on the classical sheet (see 
Figure 1). Using (3.4), we see that the function ^/ has only one square root branch cut 
[oj , d\\ on the classical sheet. The decomposition (3.31) yields a similar decomposition of 
the characteristic function, 

d I 

^{z) = n =^/(^)^^ n - ^Q^" ' (3-33) 

1=1 Q=l 

with ^ 

^^(z) = //oexp / ^i{z')dz'. (3.34) 

Equation (3.33) is the quantum version of the classical formula 

d e 

1=1 Q=l 

The basic property of the expansion (3.31) is that each factor depends only on the 
eigenvalues Xi for ^ Nj + 1 < i < ^ j<j Nj. The same property is true for tlic that 
are expressed by (3.34) in terms of only. More precisely, in the quantum theory, the 
classical eigenvalue Xi is replaced by the full set of correlators {(x^), k > 1}, or equivalently 
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by a probability distribution pi. This probability distribution is the same for all i satisfying 
the condition Yl,,j<i + 1 ^ ^ ^ S j<j" ^'^'^ ^i^^ t)e denoted by p/. The resolvent 
depends only on pi and can be expressed as 

^i{z) = ^{z; pi) = f Bli^^ . (3.36) 
J z z 

Conversely, the distribution pj is given by the discontinuity across the classical sheet branch 
cut, 

Pi(x) = ^{''%i{x + ie-) -^i{x-ie)) , xe[aj,af]. (3.37) 

A potentially confusing point^ is that the expansions (2.38) and (3.31) of =^ are similar- 
looking (at least for Nf = 0), and that the forms Mjdz and ^/dz both satisfy 

ijidz = Sij, ^idz = Sij. (3.38) 

These similarities are very misleading, and the expansions (2.38) and (3.31) are of com- 
pletely different origins. In particular, the ^/ defined by (3.32) are not single-valued on 
the curve (2.31) (even though the sum (3.31) obviously is), they have only one branch 
cut [aj,af] on the classical sheet (the [aj,aj] for J 7^ / are not branch cuts), and they 
have in general an infinite number of sheets (they arc not algebraic functions). This is 
very different from the that are single-valued on (2.31), have branch cuts on [aj,aj] 
for all J on the classical sheet, and satisfy a degree two algebraic equation. Moreover, for 
the forms ipidz, the constraint (3.38) is the non-trivial statement of Poincare duality. In 
particular, the fact that the aj-period integral of ipjdz is zero for J 7^ 7 is a non-trivial 
statement (it gives constraints on the polynomials pi in (2.36)). On the other hand, the 
similar statement on ^jdz is trivial. It just follows from the fact that has no branch 
cut or other singularity on [aj, a'j] or inside aj, and thus the contour can be deformed to 
a point when doing the corresponding integral. 

The functions or are important for our purposes because they depend only on 
pi, and thus the action of the permutation of the eigenvalues can be easily studied on them 
as we shall see below. This simple property is not shared by the i/jj. 

3.4.2 Permuting the eigenvalues 

In the sake of simplicity, let us first consider the case where all the Nj = 1 and vq = 
(this is relevant, for example, to the solution of the M = 2 theories, see Section 5). The 
generalized Konishi anomalies imply that the distributions pi are parametrized by the 
branching points of the curve (2.31) and the parameters bq (the latter enters the 
solution through the disk resolvent in (2.18)), 

J^i{z) = ^{z; pi) = ^i{z; a^ibg) . (3.39) 

The permutation operator ajj = ajj exchanges the I^^ and J*'^' eigenvalues, that is to say 
exchanges the distributions pj and pj. In explicit formulas, it is straightforward to see 

thank the referee for raising this point. 
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from the solution described in Section 2 that this amounts to exchanging the branching 
points af and Oj, because we have 



pi{x; ...,af,...,aj,...;bQ) = {ajj ■ pj){x; . . . ,af , . . . ,aj , . . . ;bQ) 

= pj{x;...,aj,...,af,...;bQ). 

In particular, this shows that 

^i{z; ...,af,...,aj,...;bQ) = {ajj ■ ^j){z; . . . ,af , . . . ,aj , . . . ■,bQ) 



(3.40) 



^j{z;...,aj,...,af,...;bQ) 



(3.41) 



and 



^Kiz; ... ,aj ,. . . ,aj,. . .-bg) = (aij ■ ,^k){z; . . . , , . . . , , . . . ; 6q) 

= ^k{z; . . . , a^, . . . , , . . . ; 6q) for ivT {/, J} . 



(3.42) 



These equations are of course valid on all sheets. If ^^ifiiz) is the value of ,^i{z) on the 
classical sheet, and ,^i^i{z) the value on the second sheet that is glued to the classical sheet 
by the branch cut [aj, af] {^i has in general additional sheets, but they are irrelevant for 
our purposes), then 

.^ifl{z] . . . , , . . . , Oj , . . . ; 6q) = ^.jfl{z- . . . , a^, . . . , , . . . ; 6q) , (3.43) 
^i,i{z; ...,af,...,aj,...;bQ) = ^j,i{z; . . . , a^, . . . , af , . . . ; 6q) . (3.44) 

We can now implement the constraint (3.9) on It is equivalent to the identity 

^(z;...,af,...,af,...;6Q) = {aij ■ ^){z; . . . ,af, . . . ,a% . . . - bg) 

= ^{z;...,aj,...,aj,...;bQ), 

which must be valid on all sheets.^ On the classical sheet, this is of course automatically 
satisfied by using (3.43), 

JV 

^o{z;...,af,...,a%...;bQ) = J| ^k,o{z; . . . ,af , . . . ,a% . . . •,bQ) 

K=l [o.W) 

= ^o{z;...,aj,...,af,...;bQ). 

Let us now consider the analytic continuation ^^^"^ of through the 7*^ branch cut (a 
priori, and ^^'^^ could be distinct functions if / 7^ J). Remembering that [aj,af] is 
not a branch cut of the ^k,o for K ^ I, we find 



ril), 



K^I 



^Let us stress that the action of aij does not mix up the different sheets. The a**" sheet can be 

characterized by the coefficients Ca,k{pi) of the large z expansion, which have a certain sheet-dependent 
functional dependence on the distribution pi. Permuting pi and pj clearly does not change a. 
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The identity (3.45) then yields 



^i,i{z; . . . ,af , . . . ,aj , . . . ■,bQ) ^k,o{z; ■ ■ ■ ,af , . . . ,aj , . . . ■,bQ) = 

^i,i{z- . . . , a^, . . . , a]^, . . . ; 6q) JJ ^k,q{z; . . . ,a% . . . ,af, . . .•,bQ) . (3.48) 

K^I 

Using (3.43) and (3.44) on the right hand side of this equation, we obtain 

^i,i{z; ...,af,...,aj,...; bQ)^jfi{z; . . . , , . . . , Oj , . . . ; 6q) = 

^j,i{z; ...,af,...,aj,...; bQ),^jfi{z; . . . ,af , . . . ,aj , . . . ■,bQ) . (3.49) 

This identity simply tells us that the analytic continuation of J^q through the J^^ branch 
cut, 

^i'\z) = ^j,i{z) n '^kA^) , (3.50) 

must be equal to the analytic continuation (3.47) through the P^^ branch cut. Using (3.26), 
this immediately yields the quantization condition (3.5). 

The most general case or arbitrary Nj and i^q can be treated along similar lines. 
One must take into account that all the eigenvalues that arc sitting at the same point 
classically are constrained by the generalized Konishi equations to have the same resolvent 
(3.32) quantum mechanically. Moreover, the eigenvalues at bq cannot fluctuate. The 
constraints on ^ are then obtained by exchanging the Nj eigenvalues at aj with the Nj 
eigenvalues at aj. This is a symmetry as long as one exchanges Nj and Nj at the same 
time, corresponding to a combination of permutations that exchange the and J*^ blocks 
of the matrix X. 



4. One more Ward identity 

Taking into account (2.47), the quantization conditions (3.5) are equivalent to 

OWby _ 5Wdv 
dsj dsj 



(4.1) 



We can use these d — 1 equations to express the d parameters sj in terms of a single 
parameter s = sj, 

SI = si{s) . (4.2) 

This parameter s is the last indeterminate in our problem. The equation (2.33) shows that 
it is identified with the glucball supcrficld (2.12) expectation value, s = {S). It is natural 
to introduce a new analytic function defined by 

WMs) = Wby{si{s)). (4.3) 

To complete the proof of the Dijkgraaf-Vafa conjecture (1.1), we need to show that 

Kw{s = {S))=0, (4.4) 
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where the ' denotes the derivative with respect to s. Let us note that when d = 1 and 
the low energy gauge group has a single U(A^i) factor, Wbv = W^y and equation (4.4) is 
all we have. This simple "one-cut" case contains already a great deal of interesting non- 
perturbative physics [22, 9]. The aim of this Section is to show that (4.4) follows from a 
simple Ward identity. 

4.1 The quantum glueball superpotential 

The simplest way to present the proof is to use the so-called quantum glueball superpo- 
tential W{s). This is defined to be the Legendre transform of VFiow; defined in (2.1), with 
respect to the gauge coupling constant. More precisely, in the theory with TVf flavors, the 
expectation value of the glueball superfield (2.12) is given by the equation (2.11) (with the 
trivial modification InA^^"''^^ — > 2i7rr when Nf = 2N). This equation can be inverted 
to find InA^^^^f as a function f of (S) = s. Indicating the dependence on the various 
couplings explicitly, with in particular (2.14), the glueball superpotential is then defined 

by 

W^(s;tK,lnA2^-^f) = Wio^{tK,f{s;tK)) + (in A^^-^f - /(s; t;^))^ . (4.5) 
By construction, W{s) is linear in lnA^^~^f, and satisfies the equations 

W'{s = {S))=0. (4.7) 

These equations are immediate consequences of the definition (4.5) and of (2.10) and (2.11). 
In particular, (4.7) implies that if we could prove that 

W = Wby , (4.8) 

then the last Dijkgraaf-Vafa equation (4.4) would automatically follow. We are thus going 
to derive (4.8) in the following. 

This derivation will necessitate a little bit of more work, but from what we already 
know we can immediately deduce that 

VF(s;tK,lnA2^-^0 = M^Dv(5;ix,lnA2^-^0 + Ms), (4.9) 

where h depends on s only (not on the couplings tx or A) and h{0) = 0. This last 
condition on h follows from the fact that, by construction, W and VFdv have the same 
classical, s — > 0, limit. Moreover, the dependence on InA^'^"'^* is through the simple 
linear term slnA^^"-^* on both sides of (4.9). This is clear for W from the definition 
(4.5), and it is also true for W^y, that depends on A only through the normalization factor 
(2.23) in the matrix integral. Actually, the A-dependence of ^ was chosen precisely for 
this matching to be valid. Finally, (4.1), (4.2) and (4.3) imply that 9f^lVDV = dtj^WDy. 
Then the equations (4.6) on the one hand, and (2.19) on the other hand, show that the 
ti^-dependence on both sides of (4.9) must be the same. 
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4.2 The last Dijkgraaf-Vafa equation 

Thus we have to prove that /i = 0, or more precisely that h must be a linear function, since 
a linear shift with a coupling-independent slope can be absorbed in the normalization =yK 
(or equivalently in a rcscaling of A on the matrix model side). 

Our strategy is as follows. Let us introduce the first order differential operator 

We are first going to show that matrix model identities imply that 

V ■ Wdv = Wdv • (4.11) 

This equation is just an identity that follows from the definition of Wbv- It strongly 
depends in particular on the s-dependence of the normalization factor yi^ (2.23). Put in 

another way, we have chosen this s-dependence precisely for the equation (4.11) to be valid. 
In particular, it fixes uniquely the Veneziano-Yankielowicz contributions to VFdv- Second, 
we are going to show that a gauge theory Ward identity implies that 

V-W = W. (4.12) 

Combining (4.11), (4.12) and (4.9), we obtain 

V-h = h = sh' , (4.13) 

whose most general solution is a linear function, as announced. 

4.2.1 Matrix model identities 

The simplest way to derive (4.11) is to start from the matrix integral and to perform the 
infinitesimal variations 

5si = esj , Ss = es , (4-14) 

on both sides of (2.13). The matrix model vacuum (2.16) is left invariant, because (2.17) 
implies that the n/ are left invariant. By taking into account the s-dependence of the 
normalization factor ^ given by (2.23), we obtain immediately 



-2^ + E ^^1^ = ^<(^F ^tree(X) + ^Qfm/'{X)Qf,)) - £^ 

= -Esk^ + -eQ'rnf{X)Q,)) - • 
^ ogk n^^ ■' " n 

k 



(4.15) 



Expanding in powers of 1/n as in (2.15), and using 

{{^Qfmf{X)Qf,))^^^ = sNf (4.16) 
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which is a direct consequence of the matrix model loop equations (see for instance [17]), 
we obtain 




S ) 



(4.18) 



(4.17) 



Relations similar to (4.17) have been considered in special cases [23]. Combining (4.17) 
and (4.18) with (2.18), we get 



and equation (4.11) then follows from the definition (4.3). 
4.2.2 Gauge theory identities 

The anomaly equations discussed in Section 2.3 are generalizations of more standard 
(anomalous) Ward identities associated with the symmetries of the gauge theory (2.3). 
There is however a useful non-anomalous U(1)r symmetry that yields a new constraint. If 
we denote by 9 the superspace coordinates, by A the gluino, and w any effective superpo- 
tential, the assignement of charges under this U(1)r is as follows: 



After Legendre transforming to W, we obtain (4.12), which completes the proof. 
4.2.3 The theory with two flavors 

The theory with two flavors can be treated along similar lines. 

On the matrix model side, if we denote collectively by tx the couplings in the tree-level 
superpotential (2.49), the planar partition function defined by (2.51) with the normalization 
factor ^ given by (2.53) satisfies 




(4.19) 



U(1)r 1 1 1 1 2 2. 



(4.20) 



This symmetry implies that 




(4.21) 




(4.22) 



This is the analogue of (4.17) for the model with fundamental flavors. If 




(4.23) 
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equation (4.22) yields 

V ■ #DV -Ns = VFdv (4.24) 

for the Dijkgraaf-Vafa superpotential defined by (2.52) and (4.3). 

On the gauge theory side, we use the anomalous U(1)r symmetry 



e \ XY tK w 

U(1)r 110 -2N 2 2. 

to show that 



(4.25) 



Wlow = E*^^ - = {T^^WuUW,Y)) - N{S) . (4.26) 

Legendre transforming, we obtain 

V -W - Ns = W . (4.27) 
Comparing this with (4.24), and using the analogue of (4.9), we conclude that W = Wdv- 
4.3 A note on the most general solutions to the constraints 

There is one implicit assumption that we have made so far.^ If iV/ / 0, we have taken for 
granted that the branch cut [aj,a'^] opens up, or equivalently that sj is not identically 
zero. This is necessary to obtain the full set of quantization conditions (3.5). However, 
another possibility is to let some of the s/ to be identically zero (let's say for 1 < I < S), 
even though the corresponding unbroken factors of the gauge group \]{Nj) are non-trivial. 
The remaining parameters sk for K > 6 + 1 are then fully determined by the quantization 
conditions f^^^ ^{z)dz G 2i7rZ for K,L > S that follow from the analysis of Section 3 
as well as from the additional constraint discussed above (for example (4.21)). Let us 
emphasize that the solutions so obtained are different from the Nj = solutions (which 
also have sj = 0), because the a/-period of ^{z) are constrained to be equal to Nj. 

What are we supposed to do with these additional solutions? One consistent possibility 
is to consider that they are unphysical and disregard them (this is consistent because they 
cannot be obtained by analytic continuation from the other solutions: they belong to a 
different phase). It is however natural to try to interpret them, since they are good solutions 
of our set of constraints. Clearly they cannot correspond to the usual vacua (2.8), and thus 
if they are physical they imply the existence of new vacua, that we denote for reasons to 
become clear in a moment 

|iVi, KS; ...■,Ns,KS; Ns+i, kg+i; . . . ; N^, ka; i^i; . . . ; ly^) . (4.28) 

For example, when the gauge group in unbroken, we would have, on top of the usual vacua 
\N,k), an additional vacuum \N, KS). In this vacuum, the expectation value of the glueball 
superfield is zero by definition, 

(iV,KS|5|iV,KS) = 0. (4.29) 



*I thank the referee for raising this point. 
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The state | A^, KS) has ah the required properties to be identified with the chirally symmetric 
vacuum postulated by Kovner and Shifman in [15]. The more general cases (4.28) are the 
immediate generalizations of the Kovner-Shifman vacua for unbroken U(A^i) x • • • x U(iVd) 
gauge group. 

The existence of the Kovner-Shifman vacua is controversial and has never been rigor- 
ously established. It is interesting to find that they naturally fit in the framework of the 
present paper, and that the most general consistent solution actually corresponds to the 
general states of the form (4.28). For this reason, and also because, to my knowledge, a 
convincing argument that rules them out has not appeared, we believe that these hypo- 
thetical new phases of A/" = 1 gauge theories deserve further investigations. Let us stress 
however that a very plausible possibility is that there exists additional constraints that 
restrict the physical vacuum states to be of the standard form (2.8). 

The original Kovner-Shifman vacua [15], if they exist, must belong to the strongly 
coupled regime of M = 1 gauge theories, and this makes their study particularly difficult. 
However, we want to point out that some of our generalized KS vacua (4.28) can be made 
arbitrarily weakly coupled by adjusting the parameters. This occurs when the unbroken 
gauge group has only abelian factors. For example, this is the case for a Coulomb branch 
with symmetry breaking pattern U(7V) — > U(l)^, which is also relevant to Af = 2 super 
Yang- Mills [7]. Thus in principle we may have KS vacua in A/" = 2 as well! The existence of 
weakly coupled KS vacua is particularly puzzling, because explicit instanton calculations 
are then reliable and seem to rule them out. There might remain the possibility that some 
subtleties of the semiclassical path integral have been overlooked. 

5. The confinement and mass gap properties 

In this last Section, we are going to apply the above results to the problem of the creation 
of a mass gap and of color confinement in = 1 super Yang-Mills theory. The strategy we 
use to prove confinement is to deform aN = 2 theory by adding a tree-level superpotential, 
as explained by Seiberg and Witten in [14]. Our contribution is to waive the two main 
non-trivial assumptions on which this famous reasoning is based. 

The first assumption is that the solution of the M = 2 theory given in [14], and 
its generalizations, are correct. This can now be derived, using the results of Sections 
2, 3 and 4, following the approach proposed in [7]. Many examples have been studied 
in the literature, for instance in [24]. However, the particularly interesting case of the 
A'^f = 2N theory, for which the Yang-Mills coupling does not run, doesn't seem to have 
been published. Moreover, even for the cases with A^f < 2A'^, which can be obtained 
from the iVf = 2N theory by integrating out the quarks, the analysis presented in the 
literature is often incomplete, because the last Dijkgraaf-Vafa equation (4.4) is usually 
not implemented. For these reasons, and also because the result is instructive, we briefly 
present the derivation for the A^f = 2A^ theory in the next subsection. 

The second and most important assumption is that the J\f = 1 theory has a mass 
gap. This is explained in Section 5.2, where we show how this fundamental property can 



- 26 - 



actually be derived by combining the solutions of both the J\f = 2 and M = 1 theories 
obtained previously. 



5.1 The H = 2 theory with A^f = 2N 

We consider the theory (2.3) with a degree d = N+1 tree-level superpotential parametrized 



as 

N 



= ePN{x) = ejl^x - ai), (5.1) 



f 

polynomials mJ of the form 



m 



7=1 



/'(X) = (X-m/),5f , (5.2) 



and we focus on the vacuum \Ni = 1;uq = Q) which has d = d = N and gauge symmetry 
breaking U(Ar) — >■ U(l)-^. The Riemann surface (2.31) takes the form 

C : = - pN-i{x) , (5.3) 

for a certain polynomial pn-i of degree degp7v-i < — 1. 

The N unknown coefficients in pn-i sue determined by the quantum equations of 
motion (2.44) and (4.4). As we have explained in details, (2.44) is equivalent to the fact 
that the characteristic function ^ is a meromorphic function on the curve (5.3). It has 
two poles of order N at infinity on the classical and second sheets, and 2Ar simple zeros at 
z = nif on the second sheet. Moreover, equation (4.4), or equivalently one of the equation 
(2.48) (with £ = 2N, Uq = 1 and A^^^^' replaced by e^*'^'^ since we are in the theory with 
zero P function) is equivalent to a simple asymptotic condition on If we note, as in 
Section 3, ^o{z) the value of ^ on the classical sheet and ^i{z) the value on the second 
sheet, then we have 

1™ ^^7 — T = e'''^ = q. (5.4) 

The above constraints can be straightforwardly solved (full details are given in [18]). The 
result is that the curve (5.3) must take the form 

N . 2N 



y' = ll{z- aif-pr,-^{z) = H^izf - f^z - ruf) , (5.5) 

7=1 ^ ^' f=l 

where pn-i and Hf^ are determined uniquely be the requirement that the second equality 
in (5.5) must hold. The characteristic function is given by 

nz) = ^{HN{z)+y{z)). (5.6) 

Prom this, we can deduce in particular the quantum chiral ring relations and the gauge 
theory resolvent ^{z). It is sometimes convenient to express the curve (5.5) in terms of 
the parameters ej defined by the relations 

N 

(TVpATf) = ^e^, l<p<N. (5.7) 

7=1 
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There is a great confusion in the literature when trying to express the Seiberg-Witten curve 
in terms of these physical parameters. However, knowing the quantum chiral ring it is a 
straightforward calculation. Introducing the polynomial 



N 



1=1 

we find 

Hn = En 



ENiz) = lliz-ei) (5.8) 

U}=iiz-mf) 



(1 - qf 



En 



where we take the polynomial part of the rational function on the right hand side of (5.9). 

The above solution does not depend on the parameter e in (5.1), and thus applies 
to the J\f = 2 theory that is obtained in the e — limit [7]. Strictly speaking, we have 
obtained the full set of chiral operator vacuum expectation values. The M = 2 prepotential 
that governs the low energy effective action can then be straightforwardly derived using 
Matone's relation and its generalizations [25]. Matone's relation itself is a consequence of 
superconformal Ward identities, as shown by Howe and West in [26]. The result is that the 
period matrix of the curve (5.5) yields the low energy U(l) coupling constants, as expected. 

The theory with ATf = 2N has a (3 function equal to zero, and thus it has long been 
conjectured to have a non-trivial S-duality invariance. Accordingly, the original curve for 
gauge group SU(2), derived in [14], was presented in terms of modular forms for a coupling 
Teff . It was soon realized, however, that even though Tcfj and the microscopic gauge coupling 
r were identical at the classical level, they could not be equal in the quantum theory 
[27]. Several curves were subsequently proposed for any gauge group SU(A^), all modular 
covariant with respect to some coupling t^s whose relation with r has remained unclear (see 
for example the discussion in [28]). Remarquably, our solution (5.5) has an unambiguous 
dependence in the UV coupling r that is not consistent with S-duality. This is in sharp 
contrast with the N = 1* theory, for which the same approach yields results in beautiful 
agreement with S-duality [1, 29]. 

5.2 The mass gap and confinement 

We start from the pure M = 2 theory, whose curve 

N 



1=1 



Y[{z - aif - Ak^^ (5.10) 



is the limit of (5.5) when the quark masses go to infinity. The limit of equation (5.9) 
implies in this case that 



TV 



(TVFXf) = ^a?, l<p<N. (5.11) 



1=1 



We can go to the SU(A^) theory and eliminate the irrelevant global massless U(l) by 
imposing the tracelessness condition Yli^i — 0- -^^^ generic values of the moduli a/, the 
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low energy theory is then a free abeUan U(l) theory. There are N special points on 
the moduli space where the curve (5.10) degenerates. The behaviour of the low energy 
U(l) couplings near these singularities shows that a dyon of magnetic charge one becomes 
massless there. The electric charge of the dyon can be any integer < A; < N— 1, depending 
on which singular point we consider. 

We now turn on the tree-level superpotential 

WtreeiX) = ^mX^ . (5.12) 

If m is infinitesimal, the resulting theory can be seen as a small deformation of the original 
J\f = 2 theory. However, the = 1 theory does not have a moduli space. The moduli of 
the M = 2 theory must be "frozen" when the superpotential (5.12) is turned on. In [14], 
Seiberg and Witten assumed that the moduli are frozen precisely at the singular points on 
moduli space. This is equivalent to assuming that the J\f =1 theory has a mass gap. Since 
this is an important point, let us present a simple proof. 

Suppose that the I\f = 1 theory obtained by deforming = 2 by adding (5.12) has a 
mass gap. Since the creation of a mass gap is a low energy phenomenon, we can analyze 
the problem in the framework of the low energy effective action. When m = 0, we have a 
free abelian theory. The only way such a theory can create a mass gap is through the usual 
Higgs mechanism. But this requires charged fields to be present in the effective action. 
This occurs only at the singular points on the M = 2 moduli space. Conversely, assume 
that by turning on (5.12) the moduli are frozen at the singular points. It is then possible to 
analyze directly the effect of Wtree in the effective action using symmetry arguments, and 
then to show that the massless dyon indeed condenses [14] . Since the field that condenses 
is magnetically charged, this also implies (oblique) confinement. Using analyticity, one 
can actually compute the dyon condensate for any value of m [30, 9] and show that it is 
non-zero including in the limit m ^ oo where we flow to the pure = 1 gauge theory. 

We are now reaching the finishing line. By using the solution of the M = I theory with 
the tree-level superpotential (5.12) described in Sections 2-4, we can compute exactly the 
correlators {TipX^). Using (5.11), it is straightforward to check that, indeed, these values 
correspond exactly to the singular points on the J\f = 2 moduli space. Prom what we have 
said before, this demonstrates that a mass gap is created and that confinement occurs in 
the J\f = I super Yang-Mills theory. 
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